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WEAK VERSUS P-SOLUTIONS TO LINEAR HYPERBOLIC 
FIRST ORDER SYSTEMS WITH CONSTANT COEFFICIENTS 

NIKOS KATZOURAKIS 


Abstract. We establish a consistency result by comparing two independent 
notions of generalised solutions to a large class of linear hyperbolic first order 
PDE systems with constant coefficients, showing that they eventually coincide. 
The first is the usual notion of weak solutions defined via duality. The second is 
the new notion of P-solutions introduced in the recent paper [K8], which arose 
in connection to vectorial Calculus of Variations in L°° and fully nonlinear 
elliptic systems. This new approach is a duality-free alternative to distributions 
and is based on the probabilistic representation of limits of difference quotients. 


1. Introduction 


Let n, N G N, T > 0 and consider the following archetypal system of first order 
evolutionary PDE with constant coefficients: 

(1.1) D t u + A:Bu = /, in (0,T) x R n , 

where u, / : (0, T) x R n —» R-^. Here A : R iVxn —» R^ is the linear mapping 
given in index form for any Q G R iVxn by 

N n 

(1.2) A:Q = Y, E ( A c«Q«)e Q 

a,p=l j=l 

In (1.1)-(1.2), we use the symbolisation Du(t,x) for the N x n spatial gradient 
matrix (D i M Q ,(t, x ))with respect to x G R” (where D, = d/dxi) , I), denotes 
the temporal derivative with respect to t G (0, T), {e * 1 , e N } is the standard basis 
of and we will assume that the putative solution u = (ui, un) t and the right 
hand side / = (/i,/ at) t lie in L 2 ((0,T)x R", R^). In indices, (1.1) reads 

N n 

V t u a + EE I-h^/3 — /"a, Oi — 1,...,1V. 

/3=1 i= 1 


The purpose of this paper is to establish a consistency-compatibility result by com¬ 
paring two fundamentally different notions of generalised solutions to a large class 
of systems as in (1.1). By assuming a certain hyperbolicity condition, we prove 
that both notions eventually agree for (1.1). The first one in the usual notion of 
weak/distributional solutions defined through duality, that is by requiring 

/ + Dc/) + f (j)\ = 0, 

J (0,T) xR n ’ 
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for all <fi £ C£°((0,T) xR n ). The second notion of solution has very recently been 
proposed by the author in [K8] and emerged in relation to the study of the (fully) 
nonlinear systems arising in vectorial Calculus of Variations in L°°, as well as in 
the overlapping area of nonlinear degenerate elliptic systems (see [K8]-[K11]). For 
the sake of completeness of the exposition, at the end of the introduction we discuss 
briefly the main objects associated with these modern areas. 

Our new concept is a duality-free notion of generalised solution which applies to 
general fully nonlinear PDE systems of any order. The a priori regularity required 
for this sort of solutions is just measurability and the nonlinearities are also allows 
to be discontinuous. Since we do not need to assume that putative solutions must 
be locally integrable, the derivatives a priori may not exist not even in the sense of 
distributions. 

The starting point of our notion in not based either on duality or on integration- 
by-parts. Instead, it relies on the probabilistic representation of the limits of differ¬ 
ence quotients by using Young measures, an indispensable tool in Calculus of Vari¬ 
ations, PDE theory and general topology (see e.g. [E, FL, CFV, FG, V]). However, 
the typical use of Young measures so far has been as a convergence tool, quantifying 
the failure of weak convergence due to oscillations and/or concentrations. In par¬ 
ticular, our idea is radically different from the concept of measure-valued solutions 
[DPM] and of their descendants and siblings. 

Let us motivate the idea of this new solution concept for the particular case 
of (1.1). To this end, it will be convenient to rewrite (1.1) in a slightly different 
fashion. Let Du = [Dju | Du] : (0, T) x R n —> R iVx ( 1 + ri ) be the space-time gradient 
of a putative solution u and for convenience we set x = (x’o, x) = (t, x) £ R 1+ ”. We 
may reformulate (1.1) as 

A : Du = /, in (0, T) X M", 
i N is the augmentation of the linear map A, given by 

S a p, a, ft = 1, N- i = 0 , 

A„/ 3 i, a,/3 = 1,...,iV; i = l,...,n. 


(1.3) 

where A : K Arx l 1 + n ) 


(1.4) 


A. - 


T : ((0,T) x R”) x (R n x K JVx ( 1 + n )) 


lN 


The augmentation of A sends X = [Xo|X] £ R JVx ( 1+ri ) to the vector Xo + A:X of 
R N . Since the particular structure of (1.3) does not play any role in the foregoing 
reasoning, it is perhaps less distracting to think in the generality of the system 

(1.5) T(x,u{x),Du(x)) = 0, x £ (0, T) x R”, 
where the coefficients are given by any Caratheodory mapping 

( 1 . 6 ) 

Then (1.3) corresponds to the particular linear choice F(x, rj, X) = A : X — f(x). 
Suppose that u is a strong solution, in the sense that u is in W{^ c ((0, T ) x R”, R w ) 
and satisfies (1.5) a.e. on (0,T) x R n . Let {D 1, l }h^o denote the difference quotient 
operators. By the equivalence between weak and strong derivatives, we have 

lf(x,u(x), lim D 1 ’ ft "'u(a;)') = 0, a.e. x £ (0 ,T) xl", 

\ is—} oo / 

along infinitesimal sequences (h v )^L 1 CR \ {0}. Since T is assumed to be contin¬ 
uous with respect to the gradient variable, this is equivalent to 

lim t[x, u(x), D 1 ' ,1 *'u(Y)') = 0, a.e. x£ (0 ,T) xl". 

is—} OO \ / 
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Note that the above statement makes sense if u is merely measurable, whereas the 
latter limit may exist even if the former does not. In order to represent it, we view 
the difference quotients as a family of measure-valued maps 

6 r1 ,h u : (0,T) x R n —> ^(]R JVx ( 1+n )) 

valued in the space of probability measures over a compactification M JVx ( 1 + n ) of the 
matrix space K Arx ( 1+n ). The exact manner we compactify plays no essential role 
for the notion of solution, but addition of “infinity” is necessary due to the lack of 
any bounds for the difference quotients which may not converge in any sense. This 
makes the theory genuinely nonlinear, even for linear PDE. The aforementioned 
space is the set of Young measures (for more details see Subsection 2.2 that follows). 
Since this set of Young measures is sequentially compact when equipped with the 
appropriate weak* topology, for any infinitesimal sequence h v —» 0, there exists a 
probability-valued map Vu : (0, T) x —> ^P(]R JVx ( 1+n )) such that 

(1.7) 5pi ,h„ u ^Vu in ^((0 as v -A oo, 

along perhaps a subsequence (vk)T ■ Therefore, we arrive at the following definition 
(for more details see Section 2): 

Definition 1 (Diffuse derivatives and 25-solutions to first order systems, cf. [K8]). 
Let T be a Caratheodory mapping as in (1.6) and 

£ := {E ai | a = 1, ..., N; i = 0, 1, ..., n j 

a basis of R Arx ( 1 + n ) consisting of rank-one matrices. 

(1) The set of diffuse derivatives of a measurable map u : (0, T) x R n —»• M. N 
consists of those probability-valued maps, denoted by Vu, which arise as subse- 
quential weak* limits of its difference quotients {D 1 ’ h u}h^o as i n (1-7), where D 1,/ ' 
is taken with respect to the frame £ (Subsection 2.3). 

(2) A measurable map u : (0,T) x R" —> M. N is a 25-solution to the system (1.5) 
on (0, T) x if for any compactly supported $ £ C c (R Arx l 1+ri )), we have 

(1.8) / $(X)T(x,u(x),X)d[Du(x)\(X) = 0, a.e. x£ (0 ,T) xl", 

J E"x» 

for all diffuse derivatives. 

In general diffuse gradients may not be unique for nonsmooth maps, but they 
are compatible with weak derivatives, whilst V-solutions are readily compatible with 
strong-classical solutions. For, if u is differentiable weakly (or just in measure in 
the sense of Ambrosio-Maly [AM, K8]), then Vu is unique and Vu = Sr>u a.e. on 
(0, T ) x R n , thus recovering strong solutions directly from (1.8). Diffuse derivatives 
can be seen as measure-theoretic disintegrations whose barycentres are the distribu¬ 
tional derivatives (see [K8]). For further results relevant to I?-solutions and their 
applications, see [K9]-[K11], [KMo, AyK, KPa, CKP], [KP]-[KP3]. 

The main result herein is that weak solutions coincide with 25-solutions for (1.3) 
in the appropriate respective spaces, if A (given by (1.4)) satisfies the next hyper- 
bolicity hypothesis: 

^ f The orthogonal complement II := N(A)^ C R Arx ( 1+n ) 

I of the nullspace of A is spanned by rank-one matrices. 
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Evidently, the nullspace is given by N(A) = {X £ R Nx(i+n) | A : X = 0}. Deferring 
until Subsection 2.5 the exact meaning of (1.9), we may now state our main result. 

Theorem 2 (Equivalence of notions & partial regularity). Consider the system 
(1.3) and suppose A satisfies (1.9) and f £ L 2 ((0, T) xR", R w ). Then, a measurable 
map u : (0,T) x —> R w is a weak solution in the space L 2 ((0,T) x if 

and only if it is aV-solution (Definition 1) in the fibre space #' 1,2 ((0 ,T) xM",^) 
(see Subsection 2.4). 

Moreover, any (V- or weak) solution satisfies the property that the orthogonal 
projection of Du on the subspace II C R JVx ( 1 + n ) exists in L 2 and for any r)<g) a £ II 
we have Da(rj ■ u) £ L 2 ((0,T ) x R”). 

From the viewpoint of applications, the significance of our “linear consistency” 
result for the model system (1.1) lies in that it possibly opens up a new avenue of ex¬ 
ploration beyond degenerate elliptic systems, allowing to access the fully nonlinear 
vectorial hyperbolic realm of systems of Hamilton-Jacobi equations -as those aris¬ 
ing from non-zero sum differential games [E3] - which can not be studied with either 
duality or viscosity solution methods, but for which the framework of 25-solutions 
is applicable. 

From the intrinsic viewpoint of the theory, we provide further insights on the 
structure of the generalised objects complementing the observations in [K8], which 
suggest (additionally to the results obtained in the previously cited papers) that 15- 
solutions constitute an adaptable and proper duality-free theory for fully nonlinear 
vectorial problems. 

As mentioned earlier, the approach of 25-solutions emerged in the study of (fully) 
nonlinear degenerate elliptic systems and of higher order equations, particularly 
those arising in Calculus of Variations in L°°. The latter area is concerned with 
the study of variational problems for functionals of the form 

E x (u,0) := \\H(;u,Du)\\ Lao{oy u £ W^(R n ,R N ), O <s R", 

as well as of the associated “Euler-Lagrange” equations which describe their ex¬ 
trema. The scalar case N = 1 was pioneered in the 1960s by Aronsson ([Al]-[A7]). 
Nowadays, the scalar case is very well developed and the relevant single equations 
are studied in the context of viscosity solutions (for a pedagogical introduction see 
[K7, C]). The vectorial case N > 2 started much later in the early 2010s ([Kl]- 
[K6]). In the simplest possible case of the functional u i-A ||Du||x,=o(.), the associated 
PDE system governing its extrema is the so-called oo -Laplacian: 

Aoo u := (^Du ® Du + |Du| 2 JDuJ X (8) 1^ : D 2 u = 0, 

where [Dm] 1 := Proj( R / Dlt ))J-- The higher order case began even more recently 
([KP2, KP3, KMo, KPa]). In the exemplary case of u i-A ||D 2 u|| L oo(.), the associated 
PDE is the so-called oo -Polylaplacian and is fully nonlinear and of third order: 

A Icii := (D 2 u)® 3 : (D 3 u)® 2 = 0. 

In the papers cited above, several results regarding 25-solutions to the above equa¬ 
tions, their generalisations and other associated problems have been established. 
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2. Preliminaries, Young measures, Fibre spaces and hyperbolicity 

2.1. Basics. We begin with some basics which will be used throughout the rest of 
the paper. Firstly, for the sake of brevity, we will henceforth use the abbreviations 

R£ : = (0, T) xl", m := 1 + n. 

Our general measure theoretic and function space notation is either standard, e.g. 
as in [E2, EG] or else self-explanatory. The norms | • | appearing will always be the 
Euclidean ones, whilst the Euclidean inner products will be denoted by either 
on R'WR^ or by on matrix spaces, e.g. on R Nxni we have 

N n 

\X\ =(X:X ) 1/2 , = X,YeR N * n \ 

Oi= 1 2=0 

etc. The symbol will also denote higher order contractions as e.g. in (1.1), (1.2) 
and the exact meaning will be clear form the context. The standard bases on R" 1 , 
rA, jjivxni w ju denoted by {e ! ji}, {e a \a} and {e“ (g> e*| a,i} respectively. If 
the range of the indices is omitted (as we just did) and unless indicated otherwise, 
Greek indices will run in {1, ...,N} and Latin indices in {0,1,..., n}. We will denote 
vector subspaces of R Afxrai as well as the orthogonal projections on them by the 
same symbol. For example, the projection Proj n : M. Nxni — x R Arxni will be 
denoted by merely II. We will also employ a compactification of the vector space 
rAx«i_ p^g exac |. wa y we compactify is irrelevant from the PDE viewpoint, as 
long as the embedding space is compact and metrisable. For concreteness, we will 
utilise the Alexandroff 1-point compactification of R JVx " 1 which will be denoted by 
jgAxni _ gWxni u {oo}. Its metric distance will be the usual one which makes 
it isometric to the sphere of the same dimension (via the stereographic projection 
which identifies {oo} with the north pole of the sphere). Then, R 7Vxni becomes a 
metric vector space isometrically and densely contained in its compactification. We 
note that balls, norms and distances taken in R JVx " 1 will be the Euclidean ones. 


2.2. Young Measures. We collect for the convenience of the reader some basic 
facts about Young measures taken from [K8]. Consider the L l space of (strongly) 
measurable maps valued in the space of continuous functions over R JVxni , which 
we symbolise as idfRj, C'(R JVx " 1 )). For background material on this space we 
refer e.g. to [FL, FG, Ed, V]. The elements of this space are those Caratheodory 
functions if 1 : Rj x R iVx " 1 — x R (i.e. functions measurable in x for all X_ and 
continuous in X for a.e. x) which satisfy 


II <5 > IIl i (Kt. c '( rJVX " 1 )) :— 



max 

V AeS JVX " 1 


$Gr,A) 


di < oo. 


Its dual space is denoted by L™» (Rj, M (M. Nxni )) , where “M” stands for the 
space of Radon measures equipped with the total variation norm. The dual space 
consists of measure-valued maps Rj. 9 x i— > $(x) € Al(R JVxri1 ) which are weakly* 
measurable , that is, for any fixed Borel set U C R JVxni , the function [$(•)](£/) : 
Rj. —X R is Lebesgue measurable. The duality pairing is given by 


( 2 . 1 ) 


L™. (Rj, Al(R iVx " 1 )) xL‘( W£,C(R Nxni )) —x R, 
:= f [ $feY)d[i?(s)](Y)dx. 
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Definition (Young measures). The space of Young measures is the set of all 
weakly* measurable probability-valued maps RJ —► J 2, (R JVxni ). Hence, it can 
be identified with a subset of the unit sphere of L’fff, (RJ, Ad(R iVxri1 )): 

^(RJ,!^" 1 ) := {tf G L~.(RJ,X(I JVxni )) : r%) G ^(E Nxni ), a.e.ieR?}. 

We will equip (RJ, R Arxni ) with the induced weak* topology (which is metrisable 
and bounded sets are sequentially precompact). The next known facts about Young 
measures will be used systematically (for the proofs see e.g. [FG] and [K8]): 

(i) ^(RJ, R Arx ” 1 ) is a convex and sequentially weakly* compact set. 

(ii) All measurable mappings U : RJ —t R Wxri1 can be identified with Young 
measures and one such imbedding is given by the Dirac mass: U i—> 5jj- 

(iii) Let U v ,U°° : RJ —> R JVx ” 1 be measurable maps, Then, up to the 

passage to subsequences, we have U v —> U°° a.e. on RJ if and only if Su* 5u°° 
in ^'(RJ,I JVxni ), as 1 / ->• oo. 


2.3. Derivatives and difference quotients with respect to general matrix 
bases. Let u : RJ —► R w be any measurable map which we understand to be 
extended by zero on R 1+ ™ \ RJ. For any a G R" 1 \ {0} and any h G R \ {0}, the 
difference quotients of u along the direction a are symbolised as 


D IMx) ■■= 


u( x + ha) — u(x) 


x G 


K T- 


Let £ = {E m | a; i } be a basis of R Wxni consisting of rank-one matrices of the form 
E al = E a <8> E^ 1 . We will need to write the gradient and the difference quotients 
of u with respect to such a basis. To this end, we will use the next elementary fact 
of matrix algebra. 


Lemma 3 (Non-orthogonal expansions). Let {U al |a;i} and be two 

bases and (.,.) an inner product on the space M. Nxni . Then, there exists a unique 
fourth order tensor {C a ipj\ a,i,/3,j} such that 


C »i0j{F Pj ,X) E ai , X G R Nxni . 

a<i P,j 

Proof of Lemma 3. We begin by noting that there exists a unique set of linear 
functionals {L ai \a;i} C (R Arxrai )* such that X = '}2 ai L a i(Xf)E al , for any X G 
Rjvxm_ Tq write these functionals explicitly, apply to this expansion of X_ the 
orthogonal projection (with respect to any inner product, e.g. the Euclidean) 

n ai := 

to infer that L ai = (J\ a i{E al ))~ 1 W a i. Consider now the basis of 1-forms {F* ai \ a; ?'} 
in the dual space (R JVxni )*, where Ff ai := Similarly, there exist unique 

linear functionals {L p j\/3 m ,j} C (R^*” 1 )** for which we have the expansion L = 
"fTp j L* Sj (L) Fpj, for any L G (R Nxni )*. We may also represent L* p - via respective 
dual projections as L* pj = We finally set C ai/3j := L* pj (L ai ). □ 

Using Lemma 3 with both bases equal to £ and for the Euclidean inner product, 
we have the following expansion, in terms of directional derivatives: 

= C aipj (E pi : Du) E ai = C °Vi ( D s(«i ' «)) R ai . 

ot,i (3,j ot,i (3,j 
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Given h £ R\ {0}, we define the difference quotients of u (with respect to 8) as 


D 


iX 


pNxni 


D 


i,h. 


Cgipj 

a,i f3,j 


D 


1 ,h 


(. E 13 ■ u) 


E° 


Clearly, if Du exists weakly in L p for some p £ [1, oo), then D l h u —> Du (strongly) 
in L p as h —> 0. 


2.4. The fibre Sobolev space. The generalised solution concept of Definition 1 
is very weak and has to be coupled with an extra admissibility condition, which, 
following the “elliptic lines” of [K8] , we formulate as membership in a certain func¬ 
tional space of partially regular maps, adapted to the PDE. For A as in (1.4), let 
II be as in (1.9). The fibre space consists of maps differentiable only along certain 
rank-one directions of non-degeneracy. We begin by identifying the Sobolev space 
W 1,2 (Ry, R w ) with its isometric image W 1 ’ 2 (R%, K w ) into a product of L 2 spaces, 
via the mapping u H > ( u , Du): 

W 1 ’ 2 ^,®^) c L 2 (r? , x R Nxni ^ . 


We define the fibre space W 1,2 {\(associated with A) as the Hilbert space 

(2.2) r 1 ' 2 ^,^) := Proj i2(R?iR v xn) W^(R^RN) L 

with the natural induced norm (written for W 1,2 maps) 

IMI*' 1 ' 2 (K”) : = (|I u IIl 2 (R5.) + l|nDu||^ 2(R?) ) 

We recall that n denotes both the vector space as well as the orthogonal projection 
on it. By employing the Mazur theorem, (2.2) can be characterised as: 


yE 1 ’ 2 (R!^,R N ) 


(u, G(u)) G L 2 (RJ, l N xn) I 3 (u^f C W 1 ’ 2 (R^,R 1V ) : ) 
(ii v ,nDi(*') -(u,G(zt)) in L 2 (R^, K^x n), as v — X oo J 


We will call G(u) G L 2 (K^.,n) the fibre (space-time) gradient of u. By using 
integration by parts and the hypothesis (1.9), it can be easily seen that the mea¬ 
surable map G(u) depends only on u G L 2 (R!^,R n ') and not on the approximating 
sequence (for the proof see [K8]). Further, G(u) satisfies the “fibre derivative prop¬ 
erty” , which is 

p (g) a G n C K Arx " 1 =4> G(u) : (77 <g) a.) = Da(v ■ u), a.e. on Rj.. 


In general, the fibre spaces are strictly larger than their “non-degenerate” counter¬ 
parts and contain elements which are not even once weakly differentiable. 


2.5. On our hyperbolicity hypothesis and comparison with relevant no¬ 
tions. We how discuss briefly the precise meaning of our rank-one spanning as¬ 
sumption (1.9) and the relation to more conventional hyperbolicity notions. First 
note that, by standard linear algebra, n = N(A)- 1 coincides with the range R(A*) 
of the adjoint operator, given by 

A* : R n —x R JVx " 1 , r? V t A = J2 (VP A 0m ) e? ® eh 

/3 ,a,i 

Hence, given any (0aG R Arxni \ {0}, we have 

£ g> a G n = R(A*) ■<==> 3p £ R N : p T A = £ ® a. 
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which by (1.4), gives that II \ {0} if and only if there exists ij £ M. N such 

that r/ a = ao £„ and = aj£ a , f° r all a and all i > 1, whilst a o ^ 0. 

Conclusively, we infer that 


(2.3) (®aeh\{0} 


ao 


7 ^ 0, F.^A^ai = £ a > Va, Vi > 1. 


Recall now that (1.9) requires the existence of d linearly independent rank-one 
matrices {f 1 <g) a 1 ,..., £ d ® q d } with d < N spanning the subspace II. Consequently, 
in view of (2.3), the assumption (1.9) is equivalent to the next condition: 

The (possibly non-symmetric) (N x iV)-matrices Aj := J2 a p e“ <8) 
have a common set of d-many (possibly non-orthogonal) left eigenvectors 
{£*, ...,£ d } spanning a subspace of R N , with respective eigenvalues cr(Aj) 
= {af/ag, ...,af/ag} the components of the vectors {a 1 , ...,a d } C K rai . 

In the light of the above, our assumption (1.9) is not comparable to the standard 
hyperbolicity requirement of iV-many real distinct eigenvalues for the matrices Aj. 
In a sense, though, (1.9) can be seen as a “weak hyperbolicity” since it implies the 
existence of plane wave solutions only along certain directions ofM. N . If however 
N = d, then (1.9) is stronger. 



Remark 4. If the hypothesis (1.9) is satisfied, then the vectors j^ 1 , ...,£ d } are lin¬ 
early independent, spanning a d-dimensional subspace of R^. To see this, suppose 
for the sake of contradiction that £ d = iA p £ p for some non-zero {Ai,...,A s } 
and linearly independent ...,£ s }, s < d — 1, then by (2.3) it follows that 

A P £ p = (£ d ) T Aj= (^TApHAj = ^A p ((e) T Aj) = ]Ta p 

\ P—1 / P—1 P=1 


P= 1 


n 


which implies that 


p= i 


U* = o 


and therefore it follows that af'/ag = af/a g, showing the vectors {a 1 , a 2 , ...,a s } are 
co-linear, which contradicts the linear independence of the basis of II. 


A sufficient condition for (1.9) (or equivalently (2.4)) to hold is when the com¬ 
mutator of the matrices Ai,..., A rl vanishes: 

(2.5) [Aj,Aj-] := AjAj - AjA, = 0, i,j > 1. 

This implies that each Aj symmetric and R^ has an orthonormal basis of eigenvec¬ 
tors; (2.5) is always satisfied when min {n, N} = 1 and the fact that (2.5) implies 
(1.9) is the content of the linear-algebraic Lemma 9 at the end of the paper. How¬ 
ever, (2.5) in a sense trivialises (1.3) since orthogonality of the basis of eigenvectors 
implies that the system decouples to N independent single equations. However, the 
next example shows that even if d = N, (1.9) is a strictly weaker notion and does 
not force decoupling to independent single equations: 


Example 5. Let n = N = 2 and set A 2 = 2Ai = 2A, where 


A := 

'2 2 ' 
1 3 

and £ 1 := 

' 1 ' 
2 

, £ 2 := 

1 

-1 

, a 1 := 

1 

4 

, a 2 := 

~T~ 

1 








!_ 8 J 


L 2 j 
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Consider (1.1) with A : R 2x2 — y R 2 having as components the matrices {Ai, A 2 }. 
Then, by invoking (2.3), one easily confirms that the orthogonal complement of the 
nullspace of the augmentation A : R 2x3 — y R 2 is spanned by the rank-one matrices 
{<G (gia 1 ,^ 2 <S> a 2 }. However, (1.1) can not be decoupled, since it takes the form 

D t wi + 2 Di«i + 2D2W1 + D1M2 + 3D2W2 = /i) in R^, 

D t u 2 + 4 DiMi+4D 2 mi + 2Dim 2 + 6D 2 tt 2 = fi, in R^. 

3. Equivalence between weak and D-solutions 

In this section we establish our main result. 

Proof of Theorem 2. The proof consists of two lemmas. The idea of the proof is 
as follows: Firstly, by approximation and some partial regularity estimates , we show 
that a map u £ L 2 (RJ.,R jv ) is a weak solution to (1.1) if and only if the projection 
of the distributional gradient Du on II C R Wx " 1 is given by the fibre gradient G(u) 
(Lemma 6). Secondly, we use the machinery of D-solutions to characterise this 
partially regular map as a D-solution to (1.1) (Lemma 7, Remark 8). 

We begin with an algebraic observation. Let A be as in (1.2) and its augmenta¬ 
tion A as in (1.4). Then, if II C R Arx ” 1 is as in (1.9), we have 

(3.1) A:(ni) = A:I, 3c>0: |A:A|>c|IIA|, 
for all X £ R Nxni . 

Lemma 6. A map u : MJf —> R w in the fibre space (2.2) satisfies 

(3.2) A : G(u) = /, a.e. on RJ, 
if and only if it is a weak solution to (1.3) in L 2 (R^,R Ar ). 

Proof of Lemma 6. Suppose first that u is a weak solution to (1.3). By mollifying 
(in space-time) by convolution in the standard way (as e.g. in [E2]), for any e > 0 
there exist u e , f £ £ C°°(f£,T — e)xR n ,R Ar ) such that u £ —> u and f E —> f in 
L 2 {f5,T — <5) x R^R^) as e —y 0 for any 6 > £ > 0. Hence, in view of (3.1), 

(3.3) A: (n Du £ ) = / £ , on (d, T — S) x R n . 

Again by (3.1), (3.3) gives the estimate 

H n SM e || i 2 ((5!T _ 5)xR „ ) < G||/|| L 2 ((5)T _ 5)xRn) , 

which is uniform in e, S > 0. By the definition of the fibre space (2.2) and the 
above estimate together with the fact that u e —> u as e — y 0 in L 2 , we obtain that 
u £ W l,2 {Wf,M, N ) and in addition nDu e — y G(u) in L 2 . Thus, by passing to 
the limit in (3.3) as £ —y 0 and as S —>■ 0, we obtain that (3.2) holds, as desired. 
Conversely, suppose that (3.2) holds. Then, by (2.2) there exists an approximating 
sequence u v —y u with nDii" — y G(u), both in L 2 as v —y 00 . Hence, we have 

A: (nD/) - / = A^nDu^-GW) 

= o(!), 

asr-> 00 , in L 2 . By the above and (3.1), for any f> £ we have 

A: (v? ® D</>) + = £{— A:(nD/) +/}<(>= o(l), 
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as v —> oo. By passing to the limit, we deduce that u is a weak solution of (1.3), 
as claimed. The lemma ensues. □ 

The next result completes the proof of Theorem 2. 

Lemma 7. A map u : Ry, —> R w in the fibre space (2.2) satisfies (3.2) if and 
only if it is a V-solution to (1.3) (Definition 1) with respect to some matrix basis 
8 depending only on A. 


Proof of Lemma 7. We begin by supposing that (3.2) holds true. By the prop¬ 
erties of the fibre space (2.2), for any matrix £ (g) a G II C M Arx " 1 we have 

Da’' 1 )? • u) —> (£ <g> a) : G (u), in L 2 ( R£) as h -> 0. 


Now we invoke our hypothesis (1.9) and Subsection 2.3 to construct a basis 8 of 
t^x”! consisting of rank-one matrices and we will express the respective difference 
quotients of u with respect to 8. By (1.9), we have 

n = spanfjf? 1 ® E 1 ,..., E d (& j] C R Nxni 

for some d < N. By Remark 4, it follows that {E 1 ,..., E d } are linearly independent 
in R w . We define a basis on R 7Vxni as in Definition 1 in the following way: 

Step 1. We complete the orthogonal complement of {E 1 , ...,E d } with an orthonor¬ 
mal basis of TV — d vectors {E d+1 ,..., E N } to create a basis of R N . 

Step 2. For each of the E a, s, we complete its orthogonal hyperplane in R rai by an 
orthogonal basis {E (a)1 , ...,E {a)n } and set E (a)0 := E a . 

Step 3. We set 


f E a ®E {a)i , if a = l,...,d; i = 0, l,...,n, 

\ E a ®e\ if a = d + 1,..., N] i = 0,1,..., n. 


Then, by defining the difference quotients D 1,h u : Ry, —> R Nxni as in Subsection 
2.3 for the above basis, by Lemma 3 we have 

(3.4) nD^u—»G(u), in L 2 (RJ,II) as h —> 0. 

By (3.4) and (3.1) we obtain that 

(3.5) A : D 1,h u — > A : G(it), in L 2 (R£, R n ) as h -A 0. 

Further, for any fixed measurable set E C R^, with finite measure and any $ G 
C c (R Nxni ), by utilising (3.2), we may estimate 


(3.6) 


$(D llh u)(A:D hh u - /) 


LHE) 


— V\E\ ll < i > ||c(R JVXr *i) A : D 1,h u — A : G(u) 


A 2 (R£) 


Hence, (3.5) and (3.6) imply 


(3.7) $(D 1,/l u) (A:D 1,/i u — —» 0, in L^E, R N ) as h -> 0. 

Moreover, the Caratheodory function 


(3.8) 


V(x,X) 


*U)(A :K-f(xj) 


Xe(x ) 
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is an element of the space L 1 


^T) 


7Vx ” 1 )) (see Subsection 2.2), because 


I^ , ||l 1 (R5.,C(K JVx "i)) < 


\E\ [ max |$(X)A:X| ) 

Vxesup P ($) 1 v ’ 7 

+ Je\ { max II/IIl 2 (r?)- 

\Aesup P ($) '7 T 


Let now (/i„)5° C R\{0} be an infinitesimal sequence. Then, there is a subsequence 
h Vk —> 0 such that (1.7) holds as k —>■ oo (Subsection 2.2). By the weak*-strong 
continuity of the duality pairing (2.1), (3.7) and (1.7), we have that 


$(D 


1 ,h„ 


\l,h v 


(3.9) 


/) = / 

J E 

[ [_ 'l’(:,X)d[Du](X) 

JE J !" X "1 

$(X)(A:X-f) d[Vu](X), 


as fc —> oo. Then, (3.9) and (3.7) yield 

[_ ®(X)(A:X-f(x)) d[Vu(x)](X) = 0, a.e. x£E. 

J R NX "l ' ' 

Since E C is an arbitrary set of finite measure, is an arbitrary function in 
C c (K Arxni ) and 25u an arbitrary diffuse gradient (Definition 1), it follows that u is 
a 2?-solution on R^., as desired. 

Conversely, suppose that u is a 25-solution in the fibre space (2.2). Then, for any 
diffuse gradient 25 u and any $ £ C c (R Nxni ), it follows that 

(3.10) [_ $(X)d[0(x)](X) = 0, a.e. x € R£, 

JR«X"1 

where 0 : R ^ —> Xii oc (R Nxni ) is the weakly* measurable measure-valued map 
defined by the formula 

(3.11) (0(x),$) := [ $(X)(A:X-f(x))d[Du(x)](X), 

jRjVXnj V / 

for any $ £ C c (R Nxni ) and a.e. x £ Rj,. Evidently, for a.e. x £ R^, the 
measure 0(x) is absolutely continuous with respect to the restriction measure 
[25ii(x)]l R Arxni . From (3.11) and Definition 1 it follows that 0 = 0 a.e. on R^. 
This implies that a.e. on R^, the support of [25w(x)]l R Arx " 1 lies in the closed set 

JSfx := |xeR Arxni |A:X — f{x)\ = o|. 

Since $ has compact support in R Nxni , we infer that 

(3.12) [_ |$(X)| A :X-f(x) d[Vu{x)](X) = 0, 

for a.e. x £ R^. By considering again the function T of (3.8) and invoking (3.9) 
and (3.12), we deduce that 


(3.13) fc lim J $(D 1 ’^u(x))(A:D 1 ’ fe ^u(x)-/(x)) 


dx = 0. 













12 


NIKOS KATZOURAKIS 


We fix R > 0 and choose d> > Xf R ( 0 ), where B.r( 0) is the closed i?-ball of 
centred at the origin. Then, (3.13) gives 


pNxn-t 


lim 

k—t oo 


(3.14) 

for any R > 0. We set 


£n{ 


D i ,n jj. 

k u 


<r} 


A:D 1 ’ h "' k u(x) — f(x) 


da; = 0, 


E K := E n < x € 




s(0) ^ 0} 


and 

(3.15) 


T r (x,X) := 


(X, 

for \X\ 

\ Q{x), 

for X 


a R 

7 R 


where x H > 0(x) is a measurable selection of the set-valued mapping with closed 
non-empty values, given by 

E r 9 x i—► JSfx.nI fl (0) C R Nxn K 

The fact that O is a measurable selection of the above set-valued map, means that 

A : Q(x) = f(x) and \0(x)\ < R, a.e. x £ E R . 

Such selections exist for large enough R > 0 by Aumann’s measurable selection 
theorem (see e.g. [FL]), although in this specific case they can also be constructed 
explicitly. By using (3.15), (3.14) implies that 

lim [ A:T fl (x,D 1 ' fcl '*u(®)) -/(a) 

k-t OO J E R 

and by recalling (3.1), we rewrite this as 


dx = 0 


(3.16) lim f A:T fl fi,IID 1 ' , ‘‘'i=ii(i)] — /(a;) dx = 0. 

k—> OO J ER V / 

Hence, (3.16) implies that 

/ A:G(u)-/ < f A^f-.HD^u)-/ 

+ J AiT^^HD^ti) -A:G(u) 

< o(l) + |A| j T s (.,nD 1 - k '*u) -G(u) 
as k —> oo, and as a consequence we have 

/ A:G(u)-/ < |A| f T^f-.nD 1 ’'*-'^)-T fl (-,G(u)) 

+ |A| / T fl (-,G(u))-G(«) +o(l), 

Je r 


(3.17) 


as fc —► oo, for large i? > 0. Moreover, by assumption u lies in the fibre space (2.2). 
By invoking (3.5), the Dominated convergence theorem, the fact that \E\ < oo and 
(3.15), we may pass to the limit in (3.17) as k —> oo to obtain 

/ A:G(u)-/ < |A| [ T R (-,G(u)) - G(u) 

Je r Je r 
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for large R > 0. Finally, we may let R — > oo and recall the arbitrariness of the set 
E C M.J, and (3.15) to infer that (3.2) holds. The lemma has been established. □ 

The proof of Theorem 2 is now complete. □ 

Remark 8 (Functional representation of diffuse derivatives). In a sense, Lemma 7 
says that all diffuse gradients, when restricted to the subspace of non-degeneracies, 
have a “functional” representation inside the coefficients, given by G (u). If we 
decompose K JVxni as Ft ® n x , the restriction of any Vu to II is given by G(u), in 
the sense that Vu L II = <5g(u)j a.e. on This is a statement of partial regularity 
type for P-solutions: although not all of the diffuse gradient is a Dirac mass, certain 
restrictions of it on subspaces are concentration measures. 

We conclude the paper with a linear algebra result, perhaps of independent 
interest, in which we establish that the vector space II of (1.9) has an orthonormal 
basis consisting of rank-one matrices which can also be completed to an orthonormal 
basis of rank-one matrices spanning R Nxni , if the stronger hypothesis that the 
matrices Ai,...,A„ commute is satisfied. In this case, though, (1.1) decouples 
to N independent equations. Nonetheless, we still think that the result below is 
interesting due to its connections to the degenerate elliptic systems of [K8]. 

Lemma 9. In the setting of Subsection 2.5, if the matrices Ai,..,A n commute, 
K^x"! has an orthonormal basis of rank-one matrices such that N-many span the 
subspace II and the rest Nn-many of them span its orthogonal complement N(A): 

(3.18) n = span[{G a0 | a}} , N(A) = span [{E ai \ a\i> l}]. 

In addition, the basis {E m \ a;i} arises in the following way: there is an orthonor¬ 
mal basis {E 1 ,..., E N } ofR N and for each a an orthonormal basis {E^ a ' >l \ i} ofW 11 
such that E ai = E a ® E (a)i . 


Proof of Lemma 9. We begin by observing that from the definition of II, we get 

(3.19) n = |V G R Nxni | Y 0 ■ (—A: A') +Y:X = 0, VA e K Arxn |. 

By standard results in linear algebra ( [L] ), we obtain that the commutativity hy¬ 
pothesis of the (symmetric) matrices {Ai,...,A„} is equivalent to the requirement 
that there exists an orthonormal basis {rj 1 , ..., 77^} C R w which diagonalises all the 
matrices {Ai,..., A„} simultaneously, namely there is a common set of eigenvectors 
for perhaps different eigenvalues {c^ 1 , ...,c^ N } of A^. Thus, 

7 


whereas <r(Aj) = {c^ 1 ,..., c^ w }. We rewrite the above as 

(3.20) A-.(if®e l ) + (-c (i)a 77 a ) = 0, Va,/3;Vi> 1. 


We now define 

(3.21) N ai := 7 f 0 



7?“® ([- c «“,o,...,o ! m),...o] t ), 

(l+i)-position 


for all a and i > 1, and also 
(3.22) N a0 := 7 f 0 


rf 0 
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where c a := (c ^ a ,..., c^“) T is the a-th eigenvalue vector of {Ai,...,A„}. The 
definition of N m and (3.20) yield that A : N m = 0, for all a and i > 1. Hence, 
N al £ N(A). Moreover, by (3.19) and the fact that the TVn-many matrices {rj a ® 
e* | a,i} are an orthonormal basis of R Nxn , we have that 


Ye n <=> Y 0 • (—A: A) + Y:X = 0, 


X £ R Nxn 



Y 0 ■ ( —A: (rj a g>e*)) + Y: (rj a 0 e l ) = 0, Va; i > 1, 
Y 0 ■ (-c( i)a r) a ) + Y: (r/ a 0 e i ) = 0, Va; * > 1, 
[Fo|Y] : N ai = 0, Va; i> 1. 


Hence, H _L N(A) if and only if Y _L N al for all a and i > 1. Since N(A) = n x , this 
shows that N(A) = span[{iV a! | a; i > l}]. Moreover, the matrices N m spanning 
N(A) are linearly independent and hence exactly iVn-many. Indeed, we have 


N ai 





) = 6 a p ( C W“c^ 



It follows that for any a ^ /3, N m is orthogonal to . Moreover, for all a and 
j in {1,..., n}, (3.21) yields 

7V ai TV"- 7 c (. i ) a c U) a 


\N ai \ ' \N^\ v /l + ( c (i) Q )2^1 + ( c (i)«)2 


e (—i,+i) 


and hence for each a the set of matrices {IV OT | *} is linearly independent. Further, 
by (3.21), (3.22) we have that 


N' 


0:0 


: N 0i = (rj~ 


1 


) : {’f 


) 


(' ri a ■ rf) { [1, c (1)a ,..., • [ - c (i)/3 , 0,..., (W), ...0] } 


= 5 a p (— 


(- 


n.(i)P i Ada 


) 


(l+i)-position 


= o, 


for all a,/3 and i > 1. Moreover, by (3.22) we have 


N a0 : 




S a0 (l + c a ■ c?) 


and as a consequence the matrices {A“°|o;} form an orthogonal set of iV-many 
elements which is orthogonal to N(A). Since the dimension of the space is N + Nn, 
all the above together with (3.19), (3.21), (3.22) prove that n = span[{./V“ 0 | a}]. 

We now show that the basis {N m \a;i} can be modified in order to be made 
orthonormal and still consist of rank-one matrices. First note that the matrices 
spanning n are orthogonal and hence we only need to fix their length. Next, n^ 
can be decomposed as the following direct sum of mutually orthogonal subspaces 
N N 

U± = 0 span[{iV m |i = l,...,n}] =: 0 W Q . 

CK—1 Q=1 


Since 


W a = 0 span[ 






: i>l ] 
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by the Gram-Schmidt method, we can find an orthonormal basis of W a consisting 
of matrices of the form N al = 77 “ (g> cS a ^ 1 with ■ gS a ^ = Sij. Finally, we define 


E a0 

N a0 

= r] a <g> 

■ _/v a0 

E ai 

:= N ai = 

rf ® cS‘ 


1 

i 

\\A + |c a 2 

. c “. 


(a)i 


E a := 


£ R Nxn \ 

£ R Nxn \ 
£R N , 


g(“)° 


1 

\/l + l c “| 2 





G K 1+n , 


where a and i > 1. By the previous it follows that {E al \ a; i} is an orthonormal ba¬ 
sis of M Arx " 1 consisting or rank-one directions such that {E a0 \ a} span the subspace 
II and {E m \ a; * > 1} span its complement IF 1 . Moreover, E al = E a ® E^ 1 . □ 
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